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Abstract 

We consider the impact of oscillations on 1-200 TeV neutrinos expected from RX J1713. 7-3946. After a description 
of the nature of the source, we obtain a prediction for the neutrino fluxes, based on the intense gamma ray flux first 
seen by CANGAROO and recently measured by H.E.S.S. experiment. We study the effect of 3 flavor oscillations 
in detail and consider the impact on the muon flux induced by these high energy neutrinos, potentially observable 
by a neutrino telescope located in the Northern hemisphere. A detector in the Mediterranean with an effective 
area of 1 km 2 and unit detection efficiency should be able to see a signal of about 10 muons per year. 



Introduction 

The search for high energy cosmic neutrinos has a long 
history, tightly connected with the history of cosmic rays. 
The hope of a happy end is linked to neutrino telescopes 
presently in operation, in construction or in project. In 
the present work, we focus on a promising potential galac- 
tic source, the young supernova remnant RX J1713.7- 
3946, and calculate the expected muon signal using stan- 
dard techniques PE]. 

We begin by a description of this object (sect.[IJ, an d 
state the predicted neutrino flux in sect.|2 In sect.|3J we 
evaluate precisely the impact of three flavor oscillations 
on the flux. The effect of the absorption of the high 
energy part of the flux in the Earth and of the live-time 
of observation is considered in sect. 01 The muon signal 
in neutrino telescopes is studied in sect. [3] A summary 
of the result is offered in sect. EJ 

1 A potential neutrino source 

In this section, we recall the motivations of interest in 
supernova remnants, and a number of recent facts that 
suggest that a specific SNR, RX JI7I3.7-3946 (G347.3- 
0.5), is a promising source of high energy neutrinos. 

SNR are likely to act as accelerators of cosmic rays 
(CR). This suspicion was raised in 1934 [3] and convinc- 
ingly supported thirty year later 0] on the basis of ener- 
getics: if several percents of the energy injected by one 
supernova {£ ~ 10 51 erg) go in CR acceleration, the losses 
from the Galaxy can be compensated. Two monographs, 
appeared in 1984 and 1990, offer still very actual sum- 
maries of the astrophysics of cosmic rays pQ and the con- 
nections with particle physics 0. To a certain extent, the 
theory of acceleration of CR in SNR is still in evolution, 
but the generic expectations are stable: the cosmic ray 
flux at the SNR is expected to be a power law spectrum: 

F P = K- E- r (1) 

with index T = 2.0 — 2.4 and maximal energy E v max 
possibly as large as several PeV, as suggested by 'knee' of 
the CR spectrum seen with extensive air showers arrays. 
Due to galactic magnetic fields, we cannot trace back CR 
to their source directly, but we can reveal sources if CR 
interact with some dense target (cosmic beam-dumps). 
The CR would partially fill the dense region, their in- 
teractions would produce mesons, which would eventu- 
ally decay yielding observable gamma and neutrino radi- 
ation. The best case is a molecular cloud near to a young 

SNR 0il[7|. 



There are converging indications that this happens in 
one specific SNR visible in the Southern sky, RX J1713.7- 
3946. Let us recall the main points: (i) A strong X ray 
source has been discovered there jS], that is compatible 
with a core collapse SN exploded in A.D. 393 at a dis- 
tance of about 1 kpc [S]. (ii) The SNR is probably inter- 
acting with a molecular cloud, at about the same distance 
[TT)j . This was observed through the CO molecule; 21 cm 
hydrogen line observations corroborate these indications 
|llj . (Hi) A large portion of the X-radiation comes from 
the same region where the cloud happens to be, the col- 
umn density that produces X-ray absorption is compati- 
ble with the observed molecular cloud, and interesting de- 
tails are continuing to emerge pi IHl IHl [TBI IT6l H7| ITS ] . 
The interpretation of the X-rays as synchrotron radiation 
from 100 TeV electrons is compatible with ATCA radio 
observations, (iv) But most interestingly, the CANGA- 
ROO team did observe TeV gamma rays since several 
years ■ This suggested this source as a CR accelerator 
|2()j . Later it was claimed [2JJ that the only likely mech- 
anism to produce the bulk of gammas is the hadronic 
one (namely, 7r°'s from proton interactions). Recently, 
the H.E.S.S. experiment also reported on the observation 
of an intense source of gamma rays, with energy in the 
TeV- 10 TeV energy range |22J - This adds support to the 
overall picture, and (already with the first data) offers a 
precise determination of the photon flux. 

One could perhaps argue that none of the items dis- 
cussed above, taken alone, seems to be conclusive. Also, 
the interpretation outlined here has to face a number of 
controversial points: the distance of the object the 
compatibility with EGRET bound El (but see and 
the uniqueness of the hadronic hypothesis [2SI have been 
all questioned. Furthermore, the spectra of CANGAROO 
and H.E.S.S. do not agree well (see H3); this may be an 
indication that the systematic error for energy measure- 
ment of one or both experiments has an underestimated 
uncertainty (note however that CANGAROO has mea- 
sured photons from the N-W rim, while H.E.S.S. mea- 
sures the spectrum for the photons coming from a wider 
region). All these objections have to be seriously consid- 
ered. However, the fact remains that RX J1713-3946 is a 
very promising case for a cosmic beam dump, where the 
available observations seem to meet theoretical expecta- 
tions. In this case, the observed gamma radiation must 
be accompanied by neutrinos. In view of the interest in 
neutrino telescopes located the Northern hemisphere, this 
is a very important conclusion. Actually, we would dare 
to say that RX J1713-3946 is at present the most definite 
hope (although not necessarily 'the best') of a success- 
ful observation of cosmic neutrinos. For other possible 
sources of TeV neutrinos in the Galaxy, see ■ 
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As a matter of fact, there is already a specific calcu- 
lation of the neutrino signal from this source [2H|- We 
improve on this calculation in the following points: we 
consider deviations from the hypothesis T = 2, we include 
the effect of live-time of measurement and of neutrino ab- 
sorption, we describe the interactions at next-to-leading 
order (NLO) in QCD, and most importantly, we consider 
the occurrence of neutrino oscillations. 



2 Secondary gamma and neutrino 
radiation 

The connection between gamma and neutrino is described 
in (see in particular ref.[38] of Chapter VIII there, or 
tab.l of 1211) and in |2]. Here, we will follow this last 
reference quite closely, and describe the relation between 
secondary gamma and neutrino radiation using the for- 
mulae of cascade theory. 1 Assuming scaling, CR primaries 
and secondaries (photons, neutrinos and antineutrinos) 
have the same type of spectrum. So we take as injection 
proton spectrum a power law with spectral index T in the 
range 2 — 2.4 as in eq. ffl. and similarly for neutrinos. 
The photon spectrum from the cascade p — > ir° — ► 7 is 



_ax_ y(r) 
a„ r ' tp 



(2) 



where AX is the column density traversed by the protons 
and X p is the interaction length of CR. The effects of the 
7T° distribution (determined by strong interactions) are 
lumped into the spectrum-weighted momenta, Z plT o in 
this example. Similar expressions hold for neutrinos, as 
a sum of several (slightly more complicated) terms that 
describe the possible branches of the it and K cascades. 2 
In summary, the flux of the neutrinos of any species is 
just proportional to the photon flux: 



F u = k- F~ 



(3) 

lr The situation is simpler than for the Earth atmosphere, since 
the muons originating from the lcptonic decays of charged pions 
or kaons do not interact. The same is true for the photons from 
neutral pion and eta decays, since the molecular cloud is very thin in 
comparison with the radiation length Xq ~ 60 g/cm 2 for hydrogen 
(=the source is "gamma-transparent"). Indeed, a cloud of a few 
parsec and with at density of few times 100 particles/cm 3 has a 
column density around 0.005 g/cm 2 or smaller. 

2 For instance, since a v e comes from p — » n + — > fj,~*~ — * u e , 
we have F„ e = AX/\ P ■ Z p7r +{Y) ■ U+^+^V) ■ F p + where 
fn+u+v i s a function that can be found in Sect. 7.1 of the book 
Kl . and stands for the additional K + contribution, weighted 
with the branching ratio into leptons, BR = 0.635. Values of the 
spectrum-weighted momenta are obtained from fig. 5.5 of the same 
book; in this way, we introduce an error at the few % level. 
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Table 1: Fluxes of neutrinos and antineutrinos originat- 
ing from pion and kaon decay, relative to the photon flux. 
The latter is assumed to be a power spectrum with value 
of index indicated in the 1 st column. 



through a proportionality coefficient k that depends on 
the type of neutrino and on the spectral index. Numerical 
values are given in tabled In this manner, the neutrino 
fluxes are predicted in terms of the measured photon flux. 

H.E.S.S. data |22| in the range E = 1-10 TeV are well 
described by a power spectrum with T = 2.19±0.09±0.15, 
that is in good agreement with theoretical expectations. 
For our purposes, and since these measurements are going 
to be improved soon in the future, we shall limit ourselves 
to set r = 2.2, modeling the photon flux as follows: 



F-, 
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TcVcm 2 s 



(4) 



The corresponding neutrinos fluxes are: 



Fl = 7.3 x lO- 12 
F° = 7.4 x 10-- 
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(5) 



where the superscripts remind us that oscillations are 
not included. In this approximation, the flux of tau (anti) 
neutrinos at the source is expected to be negligible. 

An important question is which is the uncertainty on 
the neutrino/photon ratio. A primary cause is the uncer- 
tainty in the photon flux from tt°s. Beside experimental 
errors, it is possible that the gamma radiation has other, 
non-hadronic components; this will be better quantified 
with more data and when the morphology of the source 
will be understood in detail. The uncertainties in col- 
umn density AX disappear when we consider the ratio. 
Other causes of uncertainty include the one on hadronic 
interactions (=the spectrum- weighted momenta) and the 
neglected decay channels, but again this should have a 
weaker impact on the ratio. If we consider as an anal- 
ogy the predictions of the atmospheric neutrino fluxes 
|3()j . we are lead to believe that the neutrino fluxes we 
deduced should have an accuracy of 20% or better, at 
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least in the energy region of 1 — 10 TeV. Another way 
to argue for such an accuracy is to compare the results 
of our tab^with those in j^H]- Now, if one agrees that 
an accuracy of 20 % is reached, the effects of oscillations 
must be included, since as we will see they are of the order 
of 50%. 



3 Three flavor oscillations of SNR 
neutrinos 

Now, we pass to describe the effects of neutrino oscilla- 
tions. From the theoretical point of view, the situation 
is particularly simple, since the phases of oscillations are 
really very large: 
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(6) 



The conclusion is that we just need to consider averaged 
vacuum 3 oscillations |^ |22 02] EH 123- The expres- 
sion of the probability of flavor transformation is given 
in function of the mixing matrix Utj : 



(7) 



with £, I' = e, fi, t. The probabilities are the same for 
neutrinos and antineutrinos. After propagation, the neu- 
trino fluxes become: 



e,/i,r 



t 1 01 



(8) 



Adopting the standard decomposition of Uej [3H|, we 
can summarize the present experimental information as: 
6»i2 = 32.5° ± 2° (solar neutrinos and KamLAND), 23 = 
45° ± 10° (atmospheric neutrinos and K2K), 13 = 0° ± 
10° (CHOOZ), (5 CP = 0° - 360° (namely, we do not know 
the CP violating phase 8 CP , but it appears always with 
013 ). In a reasonable approximation, the symmetric ma- 
trix P (with elements Pw) is given by: 



0.6 




(9) 



An interesting question is which deviations we can 
expect. Let us assume that there are not main causes of 

3 The MSW 1361 effect does not modify the conclusion for two 
different reasons: (1) in the vicinity of the star, the matter potential 
is negligible in comparison to the vacuum term because of the small 
density in the molecular cloud; (2) inside the Earth, the converse 
happens; the matter potential is so large that any further oscillation 
is suppressed. See also I37| . 



systematic errors. Since the formal errors in the angles 
are quite small, it is useful to expand in linear approxi- 
mation in 12 , cos 2023 an d 0i3 around the central point 
012 = 32.5°, cos 20 23 = and 13 = 0°, getting: 



(10) 




where we define x = sin 20i2, y = e 2 3 + 613 and: 



612 = 2y/ X (l - X) • 5012 

e 23 = x /4 • cos 20 23 

da = y/x(l - xj/2 ■ (50i3 • cos5 c 



(11) 



Thus, the three uncertainty in P, respectively due to the 
angles 0i2, 2 3 and <5 CP (setting S8i 3 = 10°), are: 
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(12) 



From previous equation we see that the variations are 
rather small. The main effect when we are interested to 
muon signal is due to the latter two uncertainties. Com- 
bining them in quadrature we obtain the numerical ex- 
pression: 



SP^ = -<5P eM - ±5% 



(13) 



which means that the errors introduced by the uncer- 
tainties in the parameters of oscillations are negligible. 
(It means also that there is little hope to learn anything 
useful on 3 flavor oscillations). 

Thus we evaluate oscillations with mixing angles at 
central values. Using the fluxes in eq. J5j), we arrive at 
the following expectation for neutrino fluxes at Earth: 
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TeVcm 2 s 

F Va =4.3X10"" (^)- 22 T^i 

F Vc = 3.3 X 10 12 (j§y) TcVcm 2 s 



(14) 



They are the same within 20%. The expected flux of tau 
(anti) neutrinos is the same as the flux of muon (anti) 
neutrinos, which could lead to interesting signals. How- 
ever, in the following we focus just on the muon neutrino 
and antineutrino fluxes. They give rise to muons, thus 
offering a simple way to emphasize an observable signal. 
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4 Live-time and absorption in the 
Earth 

During the sidereal day which lasts 2 x r = 23' l 56 m 4 s , 
a neutrino telescope can observe a source only when the 
overwhelming background atmospheric muons is absent. 
In first approximation, this condition is met when the 
source is below the horizon. Taking the Earth's axis 
of rotation as z direction, and x axis in such a man- 
ner that the source is in the xz plane, the direction of 
the source is s = (cos 6,0, sin 6) {6 is the declination, 
6 = — 39°46' in our case) and the one of the telescope is 
t = (cos(/>cos(7rf/T),cos</>cos(7rf/T),sin</>) {(f) is the lati- 
tude). ANTARES has <j) = 42° 50' (that we adopt for nu- 
merical example), NEMO or NESTOR are more South, 
about <f) = 36°30' and <fr — 37° 33' respectively, whereas 
BAIKAL is more North <f> = 51° 50'. The cosine of the 
zenith angle cos 9z = s ■ t is thus: 

cos#z = sino"sin0 + coso"cos</>cos(7rf/r) (15) 

The origin of the time t = is the point of highest altitude 
(the apex), and conversely, the lowest altitude is reached 
when t = t. The source becomes observable after the 
time ro that satisfies cos 6 z {to) = 0. This condition can 
be satisfied if 90°- 1 6 \ > (f> > — (90° — that happens to 
be true for all detectors except BAIKAL. In other words, 
RX J1713-3946 is always observable for BAIKAL (it is 
always below the horizon, <fr — 5 > 90°), whereas for the 
other detectors, it is observable for a fraction of time 
fii v = 1 — tq/t. This can be written: 

arccos(— tan 5 tan^) 
JUv = 1 (16) 

For ANTARES this is 78 %, whereas for NEMO and 
NESTOR this is a bit less, 71 % and 72 % respectively. 

There is another effect that diminishes the number of 
observable events: High energy neutrinos are absorbed 
in the Earth before reaching the detector. This effect de- 
pends on the column density x seen by neutrinos. When 
cos Bz < 0, we have: 

x = -2.R e cos02 • p e (cos0^) (17) 

This varies with time according to eq. (jl5|l . Here, i?® = 
6.371 • 10 8 cm is the radius of the Earth, and p® (in 
gr/cm 3 ) is the average Earth density along the line of 
sight, obtained using the PREM model |39|. Now we can 
define the neutrino absorption coefficient a„ and its time 
average a„ as: 

a v {t,E) = i - e -NAx(t)*(E) 

dt a„(t,E) (18) 

a„{E) = ^r^r— 



l\Z 




Figure 1: Projection of the Earth in the xz-plane. The 
versor of the SNR (with direction "to SNR") is shown. 
Also shown the versor of the telescope at the times t = 
and t = t, when it lies in the xz-plane. 

where Na is the Avogadro number and a the effective 
cross section of neutrino absorption. The main part is 
due to CC interactions; NC interactions increase the ab- 
sorption coefficient by a small amount. 4 This conclusion 
is in agreement with what found in |41| . For antineutri- 
nos, the calculations are exactly the same. 

For the interactions (deep inelastic scattering) we ad- 
opt the recently calculated MRST2004 partons 02] and 
work at NLO in QCD L e t, us note that the mea- 

surements of HERA at s = AE p E e = 2M p E„ with E v ~ 
50 TeV give us confidence that we have an accurate de- 
scription of the interaction cross section in the relevant 
energy range. Small modifications due to nuclear medium 
(average rock nuclei in this section, water nuclei in the 
next one) are described using the simple prescriptions 

of @3]. 

5 Signal in neutrino telescopes 

Charged-current interactions of and produce muons 
and antimuons that can be observed underground. This 
is the simplest observable for a cosmic source of high- 

4 The NC cross section is ~ 1/3 of the CC one for relevant 
energies. At first, one could guess that <r = <j cc + cr nc = atot- 
But NC interactions differ from CC interactions since they do 
not absorb a neutrino; rather, they lower its energy ('neutrino re- 
generation'). We can estimate this effect by replacing a nc with 
(jp = cr nc (l — Zp), where setting E' = E/(l — y) we define 
Z r (E) a nc {E) = Jdy (1 - y)( r ~V da nc (E' ,y)/dy |I0|. For 
r 2.2, <rp is about 10 % of a cc . Using in the exponent of eq. 1181 
a = o"cc + or with r = 2.2, we conclude that a v (E) increases by 
about 5 %. 
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energy neutrinos, and it is known since long U^J. The 
reason why we prefer to concentrate on this observable 
is that the underwater detectors can achieve a very good 
angular resolution, perhaps better than one degree. This 
offers a very effective tool to reject the background of 
atmospheric neutrinos in neutrino telescopes. 

The number of muons and antimuons reaching an area 
A in a time of observation T is: 



/>OC 

tyx+P = fuv ■ A-T • / dE v F„„ {E v )x 

J Eth 

xY^{E v ,E th ){ l-a^(E v ) ) + -> I7 M ) 



(19) 



where E v is the energy of the neutrino at the point of 
interaction and E t h is the minimal muon energy that can 
be detected, and 'V M — ► 77^" stands for the contribution 
of the antineutrinos (same expression using antincutrino 
flux, cross section, and absorption coefficient). Whenever 
needed, we take as reference values: 



A = 1 km , T = 1 solar y, E x 



th 



50 GeV 



(20) 



Neutrino fluxes F and absorption coefficients a u are de- 
fined in eqs. (|5I14() and l|18[) . Finally, the probability to 
yield a muon can be calculated by the interaction cross 
sections and the muon range in water in the following 
manner: 

Y li =N A dE^-^{E v , R(E„,E th ) (21) 
JE th aE^ 

where Na is the Avogadro number; similarly for antineu- 
trinos. The muon range R{E^, E t h) can be obtained in- 
tegrating the equation: 



dR 
dE u 



1 



a + PE„ 



(22) 



where the dependence of a and j3 on E^ in water is taken 
from ref. p^. 5 On passing, we remark that occasionally 
the calculation of the the cross sections a cc and of the 
yields Y" M are done using the DIS formula at the leading 
order (LO), but using the partons calculated at NLO. 
This procedure is not consistent, and the cross sections 
and yields obtained in this way are overestimated by 10% 
at 20 TeV, and by 25% at 1 PeV. 

The parent spectrum (distribution of the events in the 
energy of neutrinos at the interaction point) is shown in 
Figure|21for five cases: 1. a fully 'idealized' case; 2.-4. the 
three cases when oscillations, absorption and live-time are 



5 We recall that in the approximation of constant coefficients, 
R{E^E th ) = l//31og[(l + B M /e)/(l + E th /e)] with 6 = a//3. In 
the energy range of interest this agrees at 10 % with the accurate 
result when a = 2.4 ■ I0~ 3 GeV/cm and = 2 X 1CT 6 cm -1 . 
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Figure 2: Distribution of muons + antimuons above 
Eth = 50 GeV due the fluxes of eqs. @) and \14\j , that 
pass through an area A = 1 km 2 in 1 yr. The higher 
curve is the 'idealized' case. The 3 middle ones show the 
impact of absorption ("abs."), live-time ("liv.") and 3 
flavor oscillations ("osc"). The last effect is universal, 
whereas the first 2 effects are estimated for <p — 42°50'. 
The lower curve includes all three effects. 

considered once at the time; 5. the case when all three 
effects are included. The typical energies are in the range 
1-200 TeV, after inclusion of absorption (that produces a 
downward shift). For E v = 50 GeV— 1 PeV, these effects 
reduce the number of events by: 



abs.: 0.81, liv.: 0.78, osc: 0.51 



(23) 



The impact of all these effects, and in particular the one 
of oscillations, are rather important. In particular, the 
number of events expected for case 1. is: 



N Ueal = 2Q1 



(24) 



This number is compatible with the 41 events found in 
(see fig. 1 there), when we consider that in [2HI the 
spectral index is assumed to be very hard, T — 2. But 
after the inclusion of oscillations, absorption and live- 
time, the decrease is much stronger: 



AW = 9 - 3 



(25) 



This is the main result of the present work. 

We present cumulative curves in Figure|31 These per- 
mit to rescale the above results if the maximal energy of 
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Figure 3: Cumulative distribution of the number of muons 
+ antimuons above Eth = 50 GeV. The curve displaced 
at high energies is the 'idealized' case (corresponds to the 
higher curve of fig\$) the other one includes the effects 
of absorption, live-time and 3 flavor oscillations (corre- 
sponds to the lower curve of fig&jj. 



the neutrino spectrum (=the cut of the power spectrum) 
is lower than the value we allowed, E umax — 1 PeV. 
For instance, if we limit ourself to what we know from 
photons, we could believe that the cut in the neutrino 
energy happens as early as at 5 TeV. This would mean 
a dramatic reduction factor of 0.31 to be applied to the 
number in eq. (|25|) . Instead, in the more realistic case 
when the proton spectrum is cut at the energy of the 
knee {E p , max — 3 PeV), one expects a cut for the neu- 
trino spectrum somewhere close to 250 TeV (and close 
to 0.5 PeV for gamma spectrum). This means that the 
number in eq. (|25() should be diminished, but only by 5 %. 

Finally, let us note that the numbers of eqs. I|24|l and 
(|25|1 apply to a detector with unit efficiency of detection. 6 



6 Summary and discussion 

The recent H.E.S.S. measurements support the view that 
RX J1713-3946 is a source of neutrinos with energies at 
TeV and above. Existing data already permit to predict 
the neutrino flux to a reasonable level of approximation. 
Future gamma-ray data should clarify the picture, and 

6 In real detectors, the efficiency is usually included in the "ef- 
fective area" , that is an increasing function of the energy. 



possibly reveal the extension of the power spectrum. 

We calculated the expected number and distribution 
of neutrino events in underwater neutrino telescopes from 
RX J1713. 7-3946. These calculations cannot be consid- 
ered definitive for a number of reasons (e.g., CR are as- 
sumed to be solely protons, a power law spectrum is as- 
sumed, 'neutrino regeneration' is treated in the simplest 
approximation, only a perfect detector is considered). 
Also, we did not attempt to estimate the background, 
though this was done purposely: we believe that it should 
be estimated during detector operation, and we are aware 
of a number of theoretical uncertainties (generally on high 
energy part of atmospheric neutrinos flux |4fc>| , and more 
specifically on the prompt contribution). 

However, we improved over the existing calculation 
of the neutrino signal from RX J1713-3946 [2Hj in several 
senses: we considered a deviation from strict equality 
r = 2, we treated the neutrino interactions at NLO, we 
estimated absorption in the Earth and live-time of data 
acquisition, and most importantly, we included 3 flavor 
oscillations. Our calculations, in particular eq. i|25|) . sug- 
gest that a detector located in the Northern hemisphere 
should have an effective area of ~ km 2 and/or a long data 
taking time in order to see RX J1713-3946 as a source of 
high-energy neutrinos. 

We thank for pleasant discussions and help V. Berezinsky, 
A. Butkevich, M. Cirelli, P. Desiati, S. Dugad, W. Ful- 
gione, P. Ghia, T. Montaruli, G. Navarra, I. Sokalsky, 
A. Strumia, R. Thorne, Y. Uchiyama. 
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